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ABSTRACT 


The  purpose  of  this  thesis  is  to  model  a  satellite  communication  system 
with  VSATs,  using  Spread  Spectrum  CDMA  methods  and  Forward  Error  Correction 
(EEC).  Walsh  codes  and  PN  sequences  are  used  to  generate  a  CDMA  system  and  EEC  is 
used  to  further  improve  the  performance.  Convolutional  and  block  coding  methods  are 
examined  and  the  results  are  obtained  for  each  different  case,  including  concatenated  use 
of  the  codes.  The  performance  of  the  system  is  given  in  terms  of  Bit  Error  Rate  (BER). 
As  observed  from  the  results,  the  performance  is  mainly  affected  by  the  number  of  users 
and  the  code  rates. 
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I.  INTRODUCTION 


A,  BACKGROUND 

Satellites  are  an  essential  part  of  our  daily  life,  and  they  have  a  very  large  usage 
ranging  from  Seareh  and  Rescue  Operations  to  Environmental  Monitoring.  The  widest 
use  of  satellites  is,  however,  in  communication  systems.  Satellites  can  cover  vast  areas  on 
the  world;  therefore,  they  are  the  nodes  where  all  links  pass  through  in  a  communications 
network.  Many  users  can  access  such  a  network  simultaneously  while  they  are  widely 
separated  geographically. 

New  concepts  in  modem  warfare  require  the  world  armies  to  employ  satellite 
communications  more  than  ever  before.  Since  the  military  requirements  are  different 
from  the  commercial  ones,  designing  and  implementing  a  military  satellite  system  is 
usually  more  complex  and  challenging.  Military  satellite  communication  systems  must 
meet  both  quantitative  and  qualitative  requirements  in  order  to  be  deployed  into  orbit. 
The  capacity  needed,  in  terms  of  the  number  of  the  subscribers  within  a  network  and  the 
throughput  or  the  data  rate  required,  is  considered  quantitative.  In  contrast,  qualitative 
requirements  are  quality  measures,  such  as  coverage  area,  protection,  survivability, 
control,  security,  durability  and  interoperability  [1]. 


B,  OBJECTIVE 

The  objective  of  this  thesis  is  to  design  a  new  satellite  ground  segment  component 
for  use  with  military  operations  and  analyze  its  performance  with  different  data  rates  and 
coding  structures.  Therefore,  this  work  is  mostly  related  with  quantitative  issues.  The 
ground  segment  is  designed  to  fit  in  the  description  of  a  Very  Small  Aperture  Terminal 
(VSAT),  which  requires  antenna  dimensions  less  than  2.4  m.  in  diameter  [6].  The 
designed  system  has  a  wide  bandwidth  and  uses  Code  Division  Multiple  Access  (CDMA) 
method. 
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C.  RELATED  WORK 

The  VS  AT  technology  is  evolving  and  more  research  is  being  done  to  further 
improve  the  VSATs.  Most  of  the  related  work  in  VSATs  is  about  networked  use  of 
VSATs  as  in  [2],  The  studies  in  [3]  and  [4]  are  about  coding  systems  of  VSATs.  Studies 
like  [5]  are  about  designing  and  implementing  terminals  for  VS  AT  systems  and  there  are 
studies  about  frequency  aspects  of  VSATs  and  Satellite  Communications  as  in  [6].  A 
similar  study  is  done  by  [21],  and  Chapter  III  of  this  thesis  was  written  jointly.  The  study 
in  [21]  mostly  deals  with  jamming  performance  of  a  VS  AT  system,  where  this  study  is 
about  VSAT  performances  for  different  coding  schemes. 


D,  ORGANIZATION  OF  THESIS 

In  Chapter  II,  we  introduce  the  Satellite  Communication  systems.  The  main 
elements  of  a  satellite  communication,  namely,  satellites,  ground  components, 
frequencies,  satellite  orbits  are  briefly  described.  A  detailed  explanation  of  VSATs  and 
their  applications  is  given. 

In  Chapter  III,  we  explain  the  properties  of  Walsh  functions,  the  extended 
orthogonality,  PN  sequences  and  Error  Correction  Coding  (ECC).  Different  coding 
architectures  like  block  coding,  convolutional  coding  and  concatenated  coding  and  their 
advantages  are  examined. 

In  Chapter  IV,  we  develop  the  uplink  model  for  the  satellite  communication  and 
find  the  Signal-to-Noise  plus  Interference  ratio  and  the  bit-error  probability  for  each 
different  coding  scheme. 

In  Chapter  V,  we  examine  the  results  of  the  simulations  for  different  systems  to 
find  the  optimum  scheme  for  various  operational  scenarios  and  requirements. 

In  Chapter  VI,  we  summarize  our  conclusions. 
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II.  SATELLITE  SYSTEMS  AND  VSATS 


A,  SATELLITE  SYSTEMS 

The  concept  of  using  satellites  dates  back  to  the  beginning  of  the  20*  century.  A 
few  years  after  the  launch  of  the  first  artificial  satellite  (Sputnik)  in  1957,  AT&T’s 
Telstar  satellite  successfully  completed  a  transatlantic  telecast,  proving  satellite 
communication  was  possible.  By  1965  a  commercial  satellite  (INTELSAT  1)  was  already 
in  orbit.  Today  more  than  2,700  artificial  satellites  are  in  space  [8]. 

1.  Basic  Satellite  System 

a.  The  Space  Segment 

The  space  segment  of  a  satellite  system  contains  the  satellite  and  the 
terrestrial  facilities  for  controlling  the  satellite.  These  facilities  include  Tracking, 
Telemetry  and  Command  (TT&C)  stations  and  a  satellite  control  center  in  which  all  the 
operations  associated  with  station  keeping  and  assessing  the  vital  functions  of  the  satellite 
are  performed. 

A  communications  satellite  has  two  main  duties: 

•  To  amplify  the  received  signal:  The  carrier  power  at  the  receiver  input  of 
the  satellite  is  in  the  order  of  100 IE  to  \nW.  The  power  that  is  transmitted  back  from 
the  satellite  is  in  the  order  of  10  IE  to  100  IE.  The  power  gain  is,  therefore,  100  to  130 
dB. 

•  To  change  the  carrier  frequency  to  avoid  interference:  Newer  satellites 
have  onboard  processing  capability,  which  demodulate  the  received  signal,  thus 
correcting  its  errors  and  sending  the  signal  back  to  earth  by  re-modulating  and  amplifying 
the  power.  This  method  results  in  lower  bit  error  probabilities. 

A  satellite  system’s  space  segment  components  can  be  listed  as  below: 

•  Power  Supply 
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Altitude  Control 


•  Station  Keeping 

•  Thermal  Control 

•  TT&C  Subsystem 

•  Transponders 

•  Antenna  Subsystem 

b.  The  Ground  Segment 

The  ground  segment  eonsists  of  all  the  earth  stations.  These  earth  stations 
ean  either  be  eonneeted  to  the  end  users  system  direetly,  as  in  VSATs,  or  by  terrestrial 
links.  A  typieal  ground  system  eonsists  of: 

•  Antenna  System 

•  Feed  System 

•  High  Power  Amplifiers  (HP As) 

•  Low  Noise  Amplifiers  (LNAs) 

•  Up/down  Converters 

•  Controlling  and  Monitoring  Equipment 

•  Ground  Communieations  Equipment  (modems,  eoders,  ete.) 

c.  Frequency  Allocation 

Satellite  systems  use  transmission  and  reeeption  of  radio  waves  to  perform 
their  tasks.  Sinee  the  frequeney  speetrum  is  a  limited  resouree,  eareful  planning  is 
required  to  prevent  interferenee.  The  frequeney  band  is  managed  by  the  International 
Telecommunications  Union.  Different  frequeneies  are  used  for  uplink  and  downlink.  The 
main  reason  for  this  frequeney  ehange  is  interferenee.  Separation  of  transmit  and  reeeive 
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frequencies  helps  reduce  interference  both  at  the  satellite  and  the  ground  receiver.  In  the 
notation  of  satellite  communication  frequency,  the  uplink  frequency  is  written  first.  So, 
6/4  GHz  C  band  satellite  communications  means  6  GHz  is  the  uplink  frequency  and  4 
GHz  is  the  downlink  frequency.  The  usual  frequencies  used  for  satellite  communications 
is  shown  in  Figure  1 . 
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Figure  1.  Frequency  Allocation  (From  [9]  ) 
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2. 


VSATs 


VSAT  stands  for  “Very  Small  Aperture  Terminal.”  VSATs  are  used  in  the  earth 
segment  and  antenna  sizes  for  a  typical  VSAT  differ  between  60  cm.  and  2.4  m.  in 
diameter.  A  VSAT  network  can  provide  one-way  or  two-way  data  communications, 
video  broadcast,  and  voice  communication.  These  services  can  be  categorized  as  point- 
to-point,  broadcast  or  interactive  networks  [6]. 

VSATs  provide  economical  private  communication  networks.  They  are  reliable 
and  their  ability  to  be  used  across  very  distant  geographical  locations  makes  them  very 
attractive  both  for  commercial  and  military  organizations.  VSATs  are  flexible  in 
architecture,  and  they  have  limited  and  known  components.  Therefore  diagnosing 
problems  and  maintaining  the  network  is  easy.  A  VSAT  network  can  work  with  different 
capacities  in  each  direction,  thus  the  nonsymmetrical  feature  can  be  an  advantage  and 
reduce  the  costs. 


Figure  2.  A  Typical  VSAT  Block  Diagram  (From  [14]  ) 


VSAT  systems  are  generally  configured  in  star  or  mesh  formations  and  include  a 
central  station  known  as  the  hub.  The  other  stations  are  located  in  different  regions, 
which  may  include  several  various  countries.  The  star  configuration  requires  connection 
of  one  VSAT  to  another  VSAT  through  the  hub  [6].  This  results  in  a  double-hop  circuit. 
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consequently  increasing  the  propagation  delay  and  requiring  twice  the  satellite  capacity. 
Also  the  hub  is  the  single  point  of  error  in  a  star  configuration,  so  for  critical  systems,  a 
physieally  diverse  hub  is  required.  A  mesh  configuration,  on  the  other  hand,  enables 
VSATs  to  connect  to  each  other  through  the  satellite  in  a  single  hop  [10]. 

The  important  considerations  for  a  VS  AT  system  vary  for  different  types  of 
installations.  However,  the  issues  listed  below  must  be  carefully  examined  in  any  case. 

•  Orbit  Visibility 

•  Frequency  Clearance 

•  Local  Weather  Conditions 

From  the  architectural  point  of  view,  the  network  size,  network  availability,  and 
redundancy  are  the  trade-off  areas  that  determine  both  the  cost  and  the  performanee  of  a 
VS  AT  system. 


3.  Satellite  Orbits 

Satellites  and  spacecrafts  that  orbit  the  earth  follow  the  same  laws  that  govern  the 
motion  of  planets  around  the  sun.  These  laws  are  known  as  Kepler’s  Laws.  They  arise 
from  Kepler’s  observation  of  movement  of  planets  around  the  sun: 

•  The  planets  move  in  a  plane;  the  orbits  described  are  ellipses  with  the  sun 
at  one  focus  (1602). 

•  The  vector  from  sun  to  the  planet  sweeps  equal  areas  in  equal  times  (the 
law  of  areas,  1605). 

•  The  ratio  of  the  square  of  the  period  T  of  revolution  of  a  planet  around  the 
sun  to  the  cube  of  the  semi-major  axis  a  of  the  ellipse  is  the  same  for  all 
planets  (1918)  [11]. 

The  third  law  of  Kepler  is  a  highly  important  one,  which  helps  us  determine  the 
geostationary  orbit  positions.  A  satellite  in  a  geostationary  orbit  has  the  same  revolution 
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period  with  earth,  and  it  always  faees  the  same  spot  on  earth.  The  third  law  can  also  be 
written  in  the  form; 


a  :  Semi-major  of  the  ellipse  formed  by  the  orbital  path 
n  :  Mean  motion  of  the  satellite  in  radians  per  second 
fi :  Earth’s  geocentric  gravitational  constant 


// =  3. 986005 xl0‘"  2 

/  sec 


(2.2) 


With  n  in  radians  per  second,  the  orbital  period  in  seconds  is  given  by: 


P  = 


In 


n 


seconds. 


(2.3) 


Therefore, 


^GSO 


'ju-p^y 

V  4;r"  , 


(2.4) 


To  find  the  distance  of  a  geocentric  orbit,  we  use  the  period  for  one  siderial  day, 
which  is  23  hours,  56  minutes,  4  seconds.  Substituting  this  period  with  the  value  in  2.2 
into  Equation  2.4,  we  find 

^G5o  “  42164  A:m. 

When  we  subtract  the  equatorial  earth  radius:  Uj,  =6378A:m,  the  geostationary 
orbital  height  is 
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^Gso  ~  ^Gso  ~^E  ~  42164 km- 637 S km  =  357S6km. 

This  value  is  generally  approximated  at  36,000  km  beeause  a  precise 
geostationary  orbit  cannot  be  attained.  The  gravitational  fields  of  the  moon  and  the  sun 
cause  a  shift  of  about  0.85°  per  year  in  inclination.  Also,  the  “equatorial  ellipticity”  of 
the  earth  causes  the  satellite  to  drift  eastward  along  the  orbit.  Therefore,  station-keeping 
maneuvers  are  performed  periodically  to  keep  the  geostationary  satellites  in  correct 
position  [12].  Geostationary  orbit  is  considered  a  natural  “resource”  and  its  use  is 
regulated  by  international  agreements. 


B,  SUMMARY 

In  this  chapter,  we  briefly  examined  the  satellite  system.  The  Ground  Segment 
and  the  Space  Segment  of  the  satellite  system  was  explained.  The  commonly  used 
frequencies  by  the  satellite  systems  were  given  and  the  mathematical  explanation  of  the 
geostationary  orbit  was  presented.  A  detailed  explanation  of  the  VS  AT  system,  which 
will  be  the  basis  of  the  further  chapters  in  this  thesis,  was  given. 
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III.  SPREAD  SEQUENCE  SYSTEMS 


Spread  Sequence  systems  are  generated  by  using  Pseudorandom  Noise  (PN) 
sequences.  Selecting  different  PN  sequences  for  different  user  groups  in  a  system 
separates  each  group  from  the  other.  These  groups  can  be  cells  in  cellular 
communications  or  VSATs  in  satellite  communication.  To  separate  the  individual  users 
in  each  group,  Walsh  functions  are  used.  When  PN  sequences  are  combined  with  Walsh 
functions,  a  system  can  be  designed  that  can  support  high  number  of  users.  When 
Forward  Error  Correction  (EEC)  is  added,  the  performance  can  further  be  increased.* 

A,  WALSH  FUNCTIONS  COMBINED  WITH  PN  SEQUENCES  [19] 

Walsh  functions  mean  everything  for  digital  DS-CDMA  communication  systems 
such  as  satellite  or  cellular  communication.  In  cellular  communication,  each  unique 
Walsh  code  in  a  Walsh  sequence  provides  orthogonal  cover  on  the  forward  traffic 
channel  within  each  cell  to  eliminate  intracell  interference.  Each  mobile  user  in  the  cell 
has  a  unique  Walsh  function  assigned,  which  encodes  the  traffic  coming  from  the  base 
station  to  the  mobile  handset.  By  applying  his  unique  Walsh  function  to  the  traffic 
coming  from  the  base  station,  the  mobile  user  despreads  and,  in  effect,  decodes  only  the 
traffic  that  is  intended  for  a  specified  user.  When  his  unique  code  is  applied  to  the 
remaining  intracell  traffic  in  the  channel,  the  orthogonality  between  user’s  code  and  other 
user  codes  zeros  out  their  interference.  A  simplified  representation  of  this  process  is 
depicted  in  Eigure  3. 


*  This  Chapter  was  written  jointly  by  the  author  and  LTJG.  Aras  [21] 
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Figure  3.  Direct  Sequence  Spread  Spectrum  (From  [19]  ) 

On  the  other  hand,  Walsh  codes  are  not  the  only  way  to  spread  the  signal. 
Pseudorandom  noise  sequences  are  used  to  spread  the  signal  as  well.  Spread 
Pseudorandom  noise  is  a  non-orthogonal  function  so  the  coverage  supplied  by  the  PN 
sequence  would  not  be  as  sufficient  as  coverage  with  a  Walsh  code.  The  intracell  traffic 
carrying  information  for  other  users  would  remain  spread  across  the  channel's  frequency 
band  causing  interference  for  the  other  users  in  adjacent  cells  due  to  imperfect 
orthogonality  between  the  cell’s  PN  sequences.  How  this  interference  affects  the 
performance  of  the  communication  system  will  be  expressed  in  detail  in  the  next  chapter. 


Figure  4.  Effect  of  PN  Sequence  on  Transmit  Spectrum 
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We  have  introduced  an  existing  example  of  Walsh  codes  and  PN  sequences 
application,  which  is  cellular  communication.  Considering  the  VSAT  system,  Walsh 
coding  and  PN  sequences  can  be  applied  to  satellite  communication  perfectly  with  some 
differences  in  modeling.  While  every  user,  in  a  cellular  communication,  communicating 
with  his  or  her  own  base  station  in  the  same  cell  with  different  Walsh  codes  and  the  same 
PN  sequence,  every  VSAT  communicates  with  a  single  communication  satellite.  In  other 
words,  in  a  VSAT  model  every  terminal  represents  cells,  which  are  separated  by  different 
PN  sequences.  Every  channel  on  the  single  VSAT  will  be  representing  users,  which  are 
separated  by  different  Walsh  Codes. 

1.  Walsh  Codes 

Walsh  functions  have  many  desirable  properties,  which  we  will  examine  in  this 
part.  We  will  look  at  autocorrelation  functions,  and  orthogonality  property,  which  may 
prove  useful  to  the  communications  engineer. 

a.  Properties  of  Walsh  Function 

Walsh  functions  are  generated  from  generator  matrices,  Hadamard 
matrices,  Rademacher  functions  and  Walsh  binary  index.  Orthogonal  Walsh  functions  are 

defined  in  order  of  N  as  ITJv  =|w^(t);tG  (0,r),7  =  0,1,...,  A-lj  ,  consisting  of  A  =  2^ 

elements  that  are  functions  of  time  and  that  have  the  following  properties  [20]. 

0  takes  on  the  values  {+1,-1}  except  at  a  finite  number  of 

points  of  discontinuity;  where  it  is  defined  to  be  zero. 

0  =  1  for  7  =  0. 

0  Wj{t)  has  precisely  j  sign  changes  in  the  interval  (0,7). 

{  f  0,  if  j 

0  \w{t)w^{t)dt  =  \  (3.1) 

0  rfj=k 
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In  Figure  5,  set  of  Walsh  functions  of  order  eight  has  been  illustrated  for 

^8  =  {%  (0,  (0,  (0,  >^3  (0,  ^4  (0,  >^5  (0,  ^6  (0,  ^7  (O} 


In  this  case,  the  interval  (0,  T)  has  been  broken  into  =  8  pieces,  each 
section  being  =  T / N  =  T /S  time  units  long. 
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Figure  5. 


The  most  significant  property  of  the  Walsh  functions  as  applied  to  the  DS- 
CDMA  communication  systems  is  that  orthogonality  between  two  different  Walsh 
sequences.  We  will  use  this  property  frequently  in  our  performance  analysis.  Details 
about  analysis  will  be  introduced  in  the  next  section.  The  orthogonality  property  implies 
that  the  integration  of  the  product  of  any  two  different  Walsh  sequences  over  a  period  is 
always  zero.  An  explanation  of  the  forth  property  is  depicted  below. 

W3(0  =  [l  1  0  0  1  1  0  0] 

w,{t)  =  [\  0  10  0  1  0  1] 

W3(0.>Vg(0  =  [0  110  1  0  0  1] 

Values  of  "-1"  in  Figure  5  indicate  logic  zero  in  this  example  where  logic 
one  is  indicated  with  "+1".  The  product  of  is  calculated  by  "exclusive  or  " 

(XOR)  function.  XOR  function  outputs  logic  one  for  the  different  inputs  and  logic  zero 
for  same  inputs.  Since  there  are  equal  numbers  of  ones  and  zeros  in  the  product  matrix  of 
,  integration  over  a  period  is  always  zero. 


b.  Extended  Orthogonality 

Orthogonality  may  not  always  be  maintained  with  three  Walsh  sequences 
from  a  set.  The  result  of  the  integration  of  two  different  Walsh  sequences  may  not  be 
orthogonal  to  the  third  Walsh  sequence. 

lo  ^  (O^j  it)Wk  (0  =  {t)w^  (0  =  T  (3 .2) 

An  example  for  the  situation  given  in  equation  (3.2)  is  given  below  for  the 
Walsh  function  set  in  Figure  5.  Notice  that  the  integration  of  over  a 

period  is  "T" 
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w,{t)  =[10100101] 
w,{t)  =[10101010] 
w,{t).w,{t)  =[0  0  0  0  1  1  1  1] 
w,{t)  =[1  1  1  1  0  0  0  0] 
w,{t).w,{t).w,{t)  =  [\  1111  1  11] 

One  way  around  this  problem  is  to  eliminate  some  of  the  Walsh  sequenees 
from  the  eomplete  set  of  funetion  and  foree  the  equation  (3.3)  to  be  valid  for  every 

Walsh  function  in  the  reduced  set  ,  which  we  use  to  define  extended  orthogonality. 

jo  w,  {t)Wj  (0  =  j^  w,  {t)wi^  (0  =  0  (3.3) 

Besides  the  advantage  of  establishing  orthogonality  for  every  channel  on  VSAT,  reducing 
the  number  of  Walsh  sequences  in  a  set  also  reduces  the  number  of  channels  [19].  In 
extended  orthogonality,  we  will  be  limiting  the  number  of  channels  per  VSAT.  In  the  real 
world  some  VSAT  terminals  may  require  more  channels  than  others  do  in  overall 
network.  In  such  cases,  a  greater  order  of  Walsh  functions  must  be  used  to  achieve 
enough  number  extended  orthogonality  Walsh  sequences.  Therefore,  the  system  can 
support  a  sufficient  number  of  channels.  The  primary  central  hubs  in  the  star 
configuration  or  the  auxiliary  hubs  definitely  need  more  complex  Walsh  functions  than  a 
remote  terminal.  Methods  for  generating  a  maximum  set  of  extended  orthogonality  Walsh 
functions  are  given  in  reference  [19]. 

c.  Autocorrelation  of  Walsh  Function 

In  this  section,  we  will  examine  the  autocorrelation  functions  for  the 
Walsh  functions  and  average  autocorrelation  for  a  set  of  Walsh  functions.  In  order  to 
explore  the  autocorrelation  function  of  Walsh  functions,  we  will  extend  our  definition 
and  consider  each  Walsh  function  to  be  periodic  with  a  period  of  T ,  which  is  consistent 
with  their  use  in  practice.  Each  channel  on  VSAT  has  a  specific  Walsh  sequence.  A 
satellite  onboard  processor  encodes  the  data  by  applying  the  channel's  entire  Walsh 
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sequence.  The  Walsh  functions  are  periodic  signals  with  a  period  of  T .  Accordingly,  the 
Walsh  function  and  the  corresponding  autocorrelation  function  is  also  periodic  with 
period  T . 

We  define  the  normalized  autocorrelation  function,  for  any 

periodic  Walsh  function  by 

1  ^ 

C(i{u)  =  -\  w,  (t)w^  {t  -  u)dt.  (3 .4) 

i  n 


Although  Equation  (3.4)  is  used  for  continuous  signals,  we  can  consider 
Walsh  sequences  as  if  they  were  discrete  signals.  This  allows  us  to  generate  the 
autocorrelation  of  the  Walsh  functions  easily.  In  spread  spectrum  systems,  the  chip 
duration,  7). ,  is  defined  to  be  the  bit  duration  divided  by  the  number  of  chips  per  bit. 


L  = 


T 

N 


(3.5) 


We  will  consider  7)  as  the  time  increments  in  discrete  Walsh  sequences  to  build  the 

autocorrelation  functions.  The  Walsh  function  presented  in  Figure  5  is  now  presented  as 
the  autocorrelation  function  in  Figure  6. 
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Figure  6.  Set  of  Autoeorrelation  Funetion  for  JVg  (From  [19]  ) 


The  superposition  of  all  the  autoeorrelation  sequenees  in  Figure  6  divided 
by  N  gives  the  "Average  Autoeorrelation  Funetion"  as: 

=  (3.6) 

^  i=0 
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where  cr.(w)was  defined  in  (3.4).  The  average  autocorrelation  function  of 


Wg  is  generated  from  (3.6). 


Figure  7.  Average  Autocorrelation  Function,  Ag(u)  for  the  Set  (From  [19]  ) 


We  find  that  the  form  of  Ag(u)  is  similar  to  the  form  of  the  autocorrelation 

function  for  a  random  binary  signal.  A  random  binary  signal  is  generated  from  an 
infinitely  long  binary  random  sequence  in  which  the  bits  are  independent  and  identically 
distributed  random  variables.  The  resulting  normalized  autocorrelation  function  for  the 
random  binary  signal  is  defined  after  [20]  by 


Aw)  = 


H 

\u\ 

otherwise 


2.  PN  Sequences 


(3.7) 


A  pseudorandom  noise  (PN)  sequence  is  a  sequence  of  binary  numbers,  which 
appears  to  be  random,  but  is  in  fact  perfectly  deterministic.  The  PN  spreading  on  the 
uplink  provides  near-orthogonality  and,  hence,  minimal  interference  between  signals 
from  each  channels  sets  from  different  VSATs.  This  allows  reuse  of  the  band  of 
frequencies  available,  a  major  advantage  of  the  CDMA.  The  Direct  Sequence  Spread 
Spectrum  (DSSS)  uses  a  secondary  modulation,  faster  than  the  information  bit  rate,  to 
spread  the  frequency  domain  content  over  a  larger  band.  The  spreading  process  not  only 
assures  the  system  to  be  more  jammer  resistant,  but  also  ensures  it  not  to  be  detected  by 

enemy  interceptors.  The  Low  Probability  Intercept  (LPI)  systems  use  spread  spectrum  as 
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well.  In  DSSS,  each  data  bit  is  modulated  by  a  Pseudo  Noise  (PN)  sequence  that 
accomplishes  the  spectral  spreading.  The  PN  sequence  consists  of  random-like  plus  and 
minus  ones,  which  are  called  “chips.”  Each  data  bit  is  modulated  with  at  least  11  to  16 
chips.  Therefore  modulated  data  bits  seem  like  PN  sequences  in  time  domain  after  being 
modulated. 


a.  Properties  of  PN  Sequences 

The  ratio  of  data  bit  intervals  to  the  chip  durations  are  known  as 
processing  gain.  The  higher  the  processing  gain,  the  better  the  autocorrelation  properties, 
and  hence,  the  better  the  ability  to  reject  narrowband  interference.  In  Figure  8-a,  a  PN 
sequence  modulates  the  data  bits,  10  10,  (Figure  8-b),  resulting  in  the  modulated  data 
bits  in  Figure  8-c. 

I  I  I  I  I  I  I  I  I 

I  >  *  *  »  I  I  *  *  *  I 

b)^  ^ _ I  _ - 

I 

I  I  i  I  I  I  I  \  p  ■ 

Figure  8.  Modulation  with  PN  Sequence 

In  this  example  the  processing  gain  is  calculated  as  13  because  there  are 
thirteen  chips  per  bit.  The  demodulation  process  is  established  in  the  satellite  receiver  by 
applying  the  exact  PN  sequence  at  the  correct  time  (Figure  9).  The  timing  process  will  be 
introduced  in  the  following  section. 
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Figure  9.  Correct  Demodulation  with  PN  Sequence 


Data  bits  (10  10)  were  recovered  after  the  applying  of  the  PN  sequence  to 
the  modulated  data  bits.  Figure  8  illustrates  the  time  domain  characteristics  of  a  DS 
Spread  Spectrum  signal.  The  frequency  domain  representation  is  given  in  Figure  10.  The 
over  modulation  of  the  data  signal  leads  to  a  lower  power  spectral  density  covering  a 
larger  frequency  band.  If  the  spread  signal  is  transmitted  in  the  presence  of  a  narrowband 
jammer,  the  despreading  operation  at  the  receiver  will  take  the  wide  band  spread 
spectrum  signal  and  collapse  it  back  to  the  original  data  bandwidth.  The  receiver  will  also 
act  on  the  narrowband  jammer  so  that  its  spectrum  is  spread  and  causes  much  lower 
interference  to  the  despread  signal.  This  is  known  as  jamming  resistance  or  the  natural 
interference  immunity  of  spread  spectrum  signals. 

We  can  see  that  DS  Spread  Spectrum  is  bandwidth  inefficient  in  that  it 
uses  N  chips  to  transmit  a  single  bit  of  information.  Without  spreading  the  spectrum,  we 
could  have  transmitted  N  bits  in  the  same  bandwidth.  This  inefficiency  is  the  tradeoff  to 
achieve  interference  rejection,  or  the  ability  to  have  reliable  communications  even  in  the 
presence  of  an  interfering  signal,  such  as  a  jammer.  It  also  reduces  the  power  spectral 
density  of  the  transmitted  signal  so  that  its  transmission  causes  less  interference  in  other 
systems  operating  at  the  same  time  on  the  same  frequency  band. 
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Recovered  Data  Signal 


Figure  10.  Direct  Sequence  Spread  Spectrum 
b.  A  utocorrelation  of  PN  Sequence 


PN  sequences  have  good  autocorrelation  properties  to  allow  the  receiver  to 
recover  bit  timing  and  to  synchronize  the  receiver  with  the  uplink.  In  satellite 
communication,  we  expect  the  users  to  transmit  spontaneously  without  any 
synchronization  between  platforms.  In  such  a  case,  a  satellite  receiver  must  apply  the 
Walsh  Sequence  and  PN  sequence  at  the  exact  time  not  to  allow  a  phase  difference 
between  the  incoming  signal  and  the  demodulation  signals.  The  satellite  receiver  decides 
the  perfect  timing  by  using  the  autocorrelation  of  the  incoming  signal.  The 
autocorrelation  function  outputs  the  peak  value  at  the  moment  of  the  correct  timing. 
Figure  8  illustrates  the  recovering  data  signal  by  applying  the  correct  PN  sequence 
without  any  phase  difference.  If  the  receiver  would  not  have  applied  the  correct  timing  to 
demodulate  the  incoming  signal,  the  data  would  not  have  been  recovered. 
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Figure  1 1 .  Incorrect  Demodulation  of  Data  Signal 
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An  illustration  of  such  a  case  is  given  in  Figure  1 1 .  In  part  (a)  is  the  PN 
sequenee  that  is  used  in  Figure  8.  Instead  of  multiplying  the  PN  sequenee  (Figure  11 -a) 
and  modulated  data  bits  (Figure  11-e),  we  multiply  the  shifted  version  of  the  PN 
sequenee  (Figure  11-b).  The  shifted  version  of  the  PN  sequenee  was  derived  by  shifting 
the  PN  sequence  one  ehip  left  (early).  The  result  was  plotted  in  Figure  1 1-d.  The  deeision 
mechanism  after  the  PN  demodulation  is  different  from  the  demodulation  with  the  Walsh 
funetions.  Due  to  perfeet  orthogonality,  the  integration  of  the  demodulated  signal  over  a 
period  has  given  either  zero  or  one  as  the  decision  statistic.  However,  the  PN  sequences 
do  not  support  full  orthogonal  coverage.  Therefore,  we  must  use  a  decision  mechanism 
other  than  the  one  in  the  Walsh  funetion.  We  eheek  the  after  demodulation-integration 
results  with  a  predefined  threshold,  which  is  zero,  to  deeide  on  the  data  bit.  Although  the 
transmitted  and  modulated  data  sequence  was  (10  10),  the  deeision  is  made  as  (0  1  0  1) 
after  the  integration  in  the  example  given  in  Figure  8. 


B,  ERROR  CORRECTION  CODING 

The  digital  communications  itself  is  prone  to  errors.  The  information  transferred 
between  the  satellite  and  the  earth  stations  is  affected  by  external  factors  and  depending 
on  the  signal  strength,  errors  are  introduced.  The  amount  of  error  that  can  be  tolerated 
depends  on  the  application.  For  voice  eommunications,  bit  error  rates  (BER)  up  to 
lO’^ean  be  acceptable,  but  data  eommunieations  require  BER  of  at  least  10”'’.  Using 
higher  signal  power  can  increase  the  performance,  but  this  itself  alone  is  not  enough  to 
correet  all  the  errors.  Also,  the  high  power  may  not  be  available  everywhere  and  the  non- 
linearities  in  amplifiers  limit  the  output  power.  Eor  these  reasons,  error  correetion 
mechanisms  are  used. 

The  principle  of  error  correetion  coding  is  adding  redundant  bits  to  the 
information  bits  and  using  these  redundant  bits  to  deteet  and  to  correet  the  errors  at  the 
reeeiver. 
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Input  Rate  (Rt,) 

Information  Bits 


Redundant  Bits 


Output  Rate  (R^,) 


Encoded  Bits 


Figure  12.  Error  Correction  Coding 


The  code  rate  is  defined  as: 

k 

r  =  — 
n 

n  :  Number  of  encoded  bits 
k :  Information  Bits 


(3.8) 


Therefore,  in  a  coded  system,  there  are  (n  -  A:)  redundant  bits.  The  coding  is 
usually  introduced  before  the  modulation  in  the  transmitter,  and  decoding  is  done  after 
the  demodulator  at  the  receiver.  In  1948,  Shannon  demonstrated  that  by  proper  encoding, 
errors  in  the  received  sequence  can  be  reduced  to  a  desired  level  without  sacrificing  the 
rate  of  information  transfer  [13].  For  an  AWGN  channel.  Shannon’s  capacity  formula  is 


(  P  ^ 

B  •  logj 

1  + 

=  B  ■  logj 

1  +  — 

1 

1  N) 

(3.9) 


C  :  Channel  capacity  (bits  per  second) 

B  :  Transmission  bandwidth  (Hz) 

P  :  Received  signal  power  (W) 

Aq  :  Single-sided  noise  power  density  (W/Hz) 


Since  P  =  ,  Equation  (3.9)  can  be  rewritten  as 


C 

B 


=  l0g2 


r 

1  + 

V 


BNj 


(3.10) 


CjB  is  the  bandwidth  efficiency. 


24 


There  are  two  methods  for  coding;  Convolutional  Coding  and  Block  Coding.  The 
block  codes  are  linear  and  calculate  the  output  frame  by  depending  on  the  current  input 
frame  only,  so  each  block  is  coded  independently  of  the  others,  and  it  has  no  memory. 
The  convolutional  codes  store  the  memory  of  previous  input  frames  and  use  this  to 
encode  the  current  input  frame.  Convolutional  codes  have  finite  memory  and  they  are 
also  linear. 


Information 

Figure  13.  A  Communication  System  Link  with  Error  Correction  Coding 

The  errors  can  be  corrected  once  they  are  detected.  The  detection  of  errors  is 
achieved  by  using  the  redundant  bits.  Detected  errors  can  be  corrected  in  several  ways. 
The  most  common  methods  are  “Error  Concealment,”  “Automatic  Repeat  Request” 
(ARQ)  and  “Eorward  Error  Correction”  (EEC). 

The  EEC  tries  to  recover  the  original  information  from  the  received  data.  A 
simple  example  of  a  (5,  2)  block  code  can  be  used  to  demonstrate  the  EEC; 


Information 

Code 

00 

00000 

01 

01011 

10 

10101 

11 

11110 

Table  1.  (5,2)  Block  Code 
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A  look-up  table  or  a  logic  construction  can  be  used  to  implement  the  values  in  the 
table.  In  Table  1,  at  least  three  bits  would  have  to  change  in  any  code  sequence  to 
produce  another  code  sequence.  Therefore,  this  code  has  a  minimum  distance  of  three. 

When  a  code  sequence  is  received,  it  is  verified  for  its  correctness.  If  there  are 
errors  in  this  received  sequence,  the  code  compares  it  with  other  sequences  and  tries  to 
correct  it.  Clearly,  the  code  can  easily  correct  1-bit  errors.  For  example,  if  the  transmitted 
sequence  was  01011  and  it  is  received  as  01001,  the  received  sequence  will  have  a 
difference  of  one  bit  from  01011,  two  bits  from  00000,  and  three  bits  from  10101  and 
11110.  Thus,  it  can  be  corrected  as  0101 1 .  On  the  other  hand,  2  bits  of  error  will  be  closer 
to  another  sequence,  and  the  decoder  will  miscorrect  the  error.  In  some  cases,  the  detector 
can  recognize  that  more  than  one  error  has  occurred  and  declare  that  there  are 
uncorrectable  errors.  The  FEC  can  be  used  with  other  correction  mechanisms  in 
conjunction,  resulting  in  even  lower  error  rates. 

The  cost  of  error  correction  coding  is  the  decrease  in  the  information  bandwidth. 
If  we  represent  the  uncoded  bit  rate  as  and  the  coded  bit  rate  as  ,  the  code  rate  is 


R, 


R 


—  =  r„ 


(3.11) 


coded 


r^is  always  less  than  1.  For  constant  carrier  power,  the  bit  energy  is  inversely 
proportional  to  bit  rate.  Therefore, 


coded 


(3.12) 


For  BPSK  modulation,  the  BER  is 


^e.uncoded  Q 


(3.13) 


Eor  a  coded  bit  stream 


P  =o 

encoded  ^ 


vv  ^0  y 


(3.14) 
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This  shows  that  is  larger  than  the  uncoded  ’  the  probability  of  bit  error 

with  coding  is  worse  than  without  coding.  However,  this  is  the  error  rate  at  the  input  of 
the  decoder.  When  the  demodulated  bit  streams  are  fed  into  the  decoder,  some  of  these 
errors  will  be  corrected  and  there  will  be  a  coding  gain. 

1.  Convolutional  Codes 

A  convolutional  code  is  generated  by  passing  the  information  bits  through  a  finite 
state  shift  register.  The  shift  register  consists  of  A  (  A:  -bit)  stages  and  n  linear  algebraic 
function  generators,  as  shown  in  Figure  14.  The  input  data  is  shifted  into  and  along  the 
register,  k  bits  at  a  time.  The  number  of  output  bits  for  each  A:  bit  input  data  sequence  is 
n  bits.  The  N  is  called  the  “constraint  length”  and  indicates  the  number  of  input  data  bits 
the  current  output  is  dependent  upon.  The  constraint  length  determines  the  complexity 
and  the  power  of  the  code  [15][16].  Convolutional  codes  can  be  described  by  their 
generator  matrices,  tree  diagrams,  trellis  diagrams  or  state  diagrams. 


-  N  stages  - ^ 


1  2 

Encoded  Sequence  (n  bits) 


Figure  14.  Convolutional  Encoder  (  From  [15]  ) 

Decoding  of  convolutional  codes  is  a  more  difficult  problem  than  encoding.  The 
function  of  a  convolutional  decoder  is  estimating  the  encoded  input  information  using  a 
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method  that  results  in  the  minimum  possible  number  of  errors.  Unlike  a  bloek  eode,  a 
convolutional  code  is  a  finite  state  machine.  Therefore,  the  output  decoder  is  a 
“maximum  likelihood  estimator”  and  optimum  decoding  is  done  by  searching  through  the 
trellis  for  the  most  probable  sequence.  Depending  on  whether  hard  decision  or  soft 
decision  decoding  is  used,  either  the  Hamming  or  Euclidian  metric  is  used,  respectively. 

Convolutional  coding  can  be  decoded  with  several  different  algorithms.  The 
Viterbi  algorithm  is  the  most  commonly  used.  Rappaport  has  given  the  following 
definition  for  the  Viterbi  algorithm,  which  is  valid  for  both  hard  and  soft  decision 
decoding  [15]; 

Let  the  trellis  node  corresponding  to  state  Sj  at  time  i  be  denoted  S. . . 

Each  node  in  the  trellis  is  to  be  assigned  a  value  V {Sj ;)  based  on  a  metric.  The 
node  values  are  computed  in  the  following  manner: 

1 .  Set  V (5*0  o)  =  0  ™d  z  =  1 

2.  At  time  i,  compute  the  partial  path  metrics  for  all  paths  entering  each 

node. 

3.  Set  V(Sj-)  equal  to  the  smallest  partial  path  metric  entering  the 
node  corresponding  to  state  Sj  at  time  i.  Ties  can  be  broken  by  previous  node 
choosing  a  path  randomly.  The  non-surviving  branches  are  deleted  from  the 
trellis.  In  this  manner,  a  group  of  minimum  paths  is  created  from  q  . 

4.  If  i  <L  +  m,  where  L  is  the  number  of  input  code  segments  (k  bits  for 
each  segment)  and  m  is  the  length  of  the  longest  shift  register  in  the  encoder,  let 
z  =  z  -I- 1  and  go  back  to  step  2. 

Once  all  the  node  values  have  been  computed,  start  at  state  time 
i  =  L  +  m,  and  follow  the  surviving  branches  backward  through  the  trellis. 

The  resultant  path  is  the  decoded  output  for  the  input  stream. 


2.  Reed-Solomon  Codes 

Reed  Solomon  (RS)  codes  are  best  known  for  their  burst  error  correction 
capabilities,  their  importance  in  concatenated  systems,  and  their  use  in  compact  disc 
audio  technology.  RS  codes  are  non-binary  BCH  codes  that  use  the  input  and  output 
alphabets,  which  have  2“  symbols.  The  block  length  of  a  RS  code  is 

2  m  1 
-1 
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RS  codes  can  correct  t  errors  within  a  block  of  n  symbols.  The  number  of  errors 
that  can  be  corrected  is  a  function  of  the  minimum  distance  of  the  code  [17]; 

'^min  =  +  1 

RS  codes  achieve  the  longest  possible  of  any  linear  code  [15]. 

RS  codes  are  decoded  by  calculating  a  syndrome  from  the  received  block  and 
known  structure  of  the  code.  The  syndromes  are  used  to  determine  an  error  locator 
polynomial,  which  is  then  solved  to  find  the  specific  error  estimates.  Once  the  errors  are 
corrected,  the  information  is  sent  out. 

3.  Concatenated  Codes 

RS  codes  are  not  very  efficient  against  randomly  distributed  errors.  However, 
together  with  convolutional  codes  in  a  concatenated  system,  they  complement  each  other 
very  nicely.  The  Viterbi  decoder  does  not  have  any  problems  accepting  soft  decisions 
from  the  channel,  and  it  delivers  short  bursts  of  errors.  The  short  error  bursts  do  not  affect 
the  RS  decoder,  as  long  as  they  are  within  a  ^bit  symbol  of  a  RS  decoder  [18]. 


D.  CONCLUSION 

This  chapter  contained  brief  explanations  of  Spread  Spectrum  Systems,  Walsh 
Sequences  and  Orthogonality,  PN  sequences  and  Error  Correction  Coding.  When  all  of 
them  are  combined,  a  robust  and  high  performance  system  can  be  designed.  The  next 
chapter  analyses  a  system  that  utilizes  all  these  components.  A  block  diagram  of  the 
system  is  given  in  Figure  15. 
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Walsh  Code  PN  Sequence 


Walsh  Code  PN  Sequence 


Figure  15.  The  Block  Diagram  of  the  overall  Satellite  Communication  System 
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IV.  UPLINK  MODEL 


A,  PERFORMANCE  ANALYSIS  OF  A  DS-CDMA  SYSTEM 

This  analysis  is  based  on  a  satellite  communication  system  with  a  total  number  of 
/VSATs  carrying  7  channels.  Although  each  VSAT  can  carry  a  certain  amount  of 
channels,  all  of  these  channels  are  not  used  at  the  same  time.  The  total  number  of  active 
channels  in  all  of  the  VSATs  cannot  exceed  the  satellite  capacity.  The  amount  of 
channels  allocated  for  each  VSAT  is  determined  according  to  the  operational 
requirements. 

The  received  information  bits  hj(t)from  a  specified  channel  (A:  =  1 )  and  a  VSAT 
( z  =  1 )  will  be  used  for  the  analysis.  In  order  to  form  the  decision  statistic,  Y,  the  received 
spread  signal  will  be  despread,  demodulated  and  integrated  over  a  one -bit  duration. 


c(f)  rrjfi)  cos{27rfj:) 

Figure  16.  Block  Diagram  of  the  Receiver 


The  received  signal  can  be  expressed  as: 

r{t)  =  S^{t)  +  n{t)  +  C{t)  (4.1) 

5*0  (t)  :  Information  +  Co-channel  Interference 
n{t)  :  Additive  White  Gaussian  Noise 
:  Inter- VS  AT  Interference 
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1.  The  Despread  Signal,  Ji(0 

The  signal  is  despread  by  applying  the  appropriate  PN  sequenee.  The  PN 
sequenee  must  be  synehronized  with  the  transmitting  VSAT.  After  dispreading,  the 
Walsh  sequenee  is  applied  to  seleet  the  desired  ehannel  and  remove  all  the  other 
ehannels’  information.  The  despread  signal,  y^{t)  is 

y,{t)  =  r{t)c{t)w,{t)  (4.2) 

Ti  (0  =  *^0  (0  c{t)  Wj  (0  +  n{t)  c{t)  (0  +  g{t)  c{t)  Wj  {t)  (4.3) 

Sq  {t)  eontains  both  information  and  eo-ehannel  interferenee: 

*^0  ~  ^  a Yd 

Ti  (0  =  A  (0  +  Yi  (0 + «i  (0 + ?i  (0 

n,  (t)  :  Additive  White  Gaussian  Noise  for  T‘  ehannel 

fj(t)  :  Inter-VSAT  Interferenee  for  P‘  ehannel 

A  it)  +  Ti  (0  =  li  ( 2PA  it)w,  it)  cos(2;rX0)  c(0w,  (0  (4.4) 

Pj^  ;  Reeeived  power  from  k*  ehannel 
h^(t)  ;  Information  bits  for  k*  ehannel  {+1,-1} 

w^(t)  :  Walsh  Code  for  k*  ehannel  {+1,-1} 

In  Equation  (4.4),  A:  =  I  is  the  index  of  the  despread  information  bits  we  want  to 
recover.  All  the  other  values  ofk{k^\)  will  be  considered  as  co-channel  interference. 

A  it)  =  ^/^  A  (0i+i  it)cit)  cos{27lfj:)cit)w,  {t) 

^  (4.5) 

=  V^A(0cos(2;rX0 
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i^-1 

S  it)Wk  (Oc(0  cos(2;rX0  c{t)w,  {t) 

)  (4.6) 

=  Z  ^k  (Ow,  {t)w^  (t)  COSilTTfj) 

k=0 

k=Al 


Additive  White  Gaussian  Noise  term: 


n,{t)  =  n{t)c{t)w,{t)  (4.7) 


Each  VS  AT  terminal  has  its  own  PN  sequence.  But  the  quasi-orthogonal  nature  of 
PN  sequences  results  in  interference.  This  interference  level  increases  as  the  number  of 
VSATs  increases.  The  VSAT  and  the  satellite  are  synchronized;  therefore,  a  phase 
difference  ( (p. ),  with  the  interfering  terminals  is  present. 


'  M  -ft-,- -I 

?i(0  =  YjH^^k  bikit +  T^i)Wikit  +  T^i)c,{t  +  ri)cos{27tfp  +  (p,)  c{t)w,{t) 


M  ft-/-! 

Z  Z  b,,  (t  +  T, )  w,,  {t  +  7-, )  w,  (0  c,  {t  +  7-,  )c{t)  cosilnfp  +  (p, ) 


2,  The  Demodulated  Signal,  yiit) 

The  demodulated  signal  is  obtained  by  applying  the  modulation  tone,  which  is 
coherent  detection. 


yiit)  =Ti(0cos(2;7X0 

=  ( A  (0  +  Ti  (0  +  «i  (0  +  ?i  (0 )  cos(2n:fj) 
=  /2(0+  T2(0  +  ^2(0  +  ?2(0 


In  order  to  obtain  the  decision  statistic,  Y,  we  need  to  calculate  each  of  the 
components  separately. 


Despread  and  modulated  information  bits,  /jCO  : 


h  =  ^ibiit)cos{27rfj)cosi27rf/) 

=  ^^^i(0(1  +  cos(4;7X0) 


(4.9) 
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where 

eos"(2;rX0  =  ^(l  +  cos(4;rX0) 

Despread  and  modulated  co-channel  interference,  5/3(0  • 


K-\ 

r2(0  =  S 


k=l 


2 


(0  (0 1^1  (0  ( 1  +  cos(4;rX0 ) 


(4.10) 


Despread  and  modulated  AWGN 

^2(0  =  n{t)c{t)w^{t)cos{27tf^t)  (4.11) 

Despread  and  demodulated  inter- VS  AT  interferenee 


M  K-\ 

?2(o=£z 


z=l  /t=0 


bik  {t  +  7-, )  (t  +  r, )  Wi  (0  c,  {t  +  r,  )c(0  cosilnf^t  +  (p, )  cos{27rfp) 


where 


1 


eos(2.;r fj:  +  (p^zo%{27t ff)  =  “[cos -  eos(4.;r fj:) eos q)^  +  sm(4.;r f/) sin ] 


M  Kf-l  I2P 

?2  (0  =  £  X  b-,  {t  +  7-. )  (t  +  T.)w,  (t)  c.  {t  +  )c{t) 


i=\  k=Q 


(4.12) 


[eos  g).  -  cos{A7r fj)  eos  +  sin(4;;r fj)  sin  ] 


M  :  Total  number  of  VSATs. 

A. :  The  total  number  of  users  in  the  /*  VS  AT. 


3.  The  Decision  Statistic 

The  deeision  statistie,  Y,  is  modeled  as  a  Gaussian  distribution  with  a  mean  value 
of  Y  and  varianee  of  ^ . 
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1 

Y  =  \y2{t)dt 
0 

T  T  T  T 

7  =  1  l2{t)dt  +  ^  Yi(f)dt  +  ^n2{t)dt  +  ^g2{t)dt 


(4.13) 


Y :  Co-channel  Interference, 


T  K-\  hp 

y-ir-  n  (0  (0(1  +  cos(4;rX0 )  dt 


(4.14) 


n  :  Additive  White  Gaussian  Noise 


1 

r]  =  \  ^(0  c(0  '^1  (0  cos(2.;r  fj)  dt 


(4.15) 


g :  Inter- VS  AT  Interference 


g  =  \g2{t)dt 

0 

M  K-\  l2p 

_  f  Z  Z  (^  +  )  ^ik  (0  c  (t + C  )c(0 

—  J  ^ 


[cos  -  cos(4;;r fj:)  cos  (p^  +  sin(4;;r fj^)  sin  (p-  ]  dt 


=  Z  Z  J  kk  (t  +  (t  +  C  )1+1  (0  c  +  C  )c(0  dt 


In  reference  [19]  a.j.(t  +  7',.)  and  (ij(t)  are  given  as 


a.k  (t  +  Tj=b„{t  +  T,  )1V.,  {t  +  T, )  C.  (t  +  r. ) 

d,{t)  =w^{t)c{t) 


(4.16) 


M  K  -X  l2p  T 

7  =  1  /f  =  n  Z  ;  ^ 


(4.17) 
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a.  Mean  Value  of  Decision  Statistic,  Y 

{r  r  r  r  r 

1 1 2  (0  dt-\-  jy2  (0  dt  +  jn2  (t)  f  2  (0  dt  +  jj2  (t)  dt 

00000 

Y  =  E\^f{t)d^  +  ^ 

=  E\^f{t)dt 

Y=  \J^bft){\  +  cos{A7ift))dt 
0  2 


2 


(4.18) 


(4.19) 


The  number  of  zeros  and  ones  are  assumed  to  be  equal,  so  the  integral  of 
bft)  over  a  period  is  equal  to  ll.Also,  f  is  a  integer  multiple  of  l/T, therefore, 
integrating  a  cosine  function  over  m  periods  is  always  zero. 


4.  Variance  of  Total  Interference  plus  Noise 


Since  the  co-channel  interference  and  noise  variances  are  independent,  the 
variance  of  total  interference  and  noise  is  the  summation  of  each  variance. 


Var[^]  =  Farjy}  +  Var\n\  +  Farjf} 


(4.20) 


a.  Variance  of  Co-channel  Interference 

Var{Y]  =  E{Y^] 

\r  _  Y' 

(0  (0  (0(1  +  cos{An  ft))  dt  > 

\S  ^  . 

=  0 


=  E\ 


Tk-x 

is 

n  k=o 


(4.21) 
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T 

Because  ^w^{t)wi^{t)dt  =  0  where  k^l. 

0 


This  is  why  we  use  Walsh  functions  for  CDMA. 
b.  Variance  of  A  WGN 

Var[n\  = 

{T  T 

j  n{t)  c{t)  Wj  {t)  cos{27I  fj)dt^  n(A)  c(A)  (A)  cos(2.;r  fA)  dA 

0  0 


{T  T 

I  j  n{t)n{A)  c{t)  c{A)w^  (t)  Wj  (A)  cos(2;;r  ft)  cos(2.;r  fA)dt  dA 

0  0 


(4.22) 


We  assume  that  n{t)  is  a  wide-sense  stationary  white  noise  process.  The 
autocorrelation  of  the  process  was  defined  in  [16]  as 


E[n{t)n{A)^  =  ^^S(t-A) . 


(4.23) 


so, 


T 

Var[n\  =  [ c^(t)wf  (t)  — (l  +  cos(4.7r ft))dt 

9  •'  9 


2 

u 

N.T 


(4.24) 


c.  Variance  of  the  Inter-  VS  A  T  Interference 

rar{g]  =  E{g^} 


M  Kf-l 


/=1  k=0 


cos  {(p.)^a,,{t  +  T.)dft)dt 


M  K--1  f  p  T  T  I 

Var{g]  =  Y,iL^\^  i^i)\  )^i  (0  dt\ a.,  (A  +  r,  )d,  (i) 


/=!  /t=0 
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M  K,-\  T  T 


Var{g]  =  ^^^\\E{a,,{t  +  T,  )a.,  (i  +  T,  )]E{d,  {t)  d,  (i)}  dtdA 

'  /=1  ^=0  0  0  ^  ^ 

'J  n/^  i\  n/^  T\ 


P(t-X) 


P(t-X) 


M  K,-\TT 


Var[g^  = 


i=\  k=0  0  0 

Applying  the  transformation  of  limits  and  variables  for  the  integral 


u  =t-A 
V  =  t  +  A 

t  =  ^{u  +  v) 
yl  =  — (v-m) 


Ja  =  det 


-^u  +  v  =  2t 

-^M-v  =  -2  A 


dt 

dA^ 

“1 

r 

du 

du 

=  det 

2 

2 

dt 

dA 

1 

1 

dv 

dv  J 

_2 

2_ 

After  transformation, 

M  ^,-1  p  T  2r-|«| 

Var{g]  =  ^'^^\  J  (d\u)J^^dudv 


1^1  ^  -T  |i,| 

M  K-\  p  t2T-\u\ 

=  I  P^(M)dudv 

i^\  4  Q  1^1 

M  Ki  -1  p  T 

i-1  fc-0  '  0 


u 


P{u)  = 


,  Nu 

' 

1 - 

\u\ <  T 

T 

0 

otherwise 
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M  K:-\  T)  /N 


/=1  k=0  “  0 


2Nu  Nu 


(2T  -2u)du 


m"  2Nu^ 


=  yy^2  Tu-Nu^+^-^  +  ^-^ 
tfh  4  [  37  2  37  47" 


^  ^(“1  p  f  27^^ 

“§§T  Uiv 


M  K.j-\  p  /  rj-i2  ' 

=  y  y  —  -  for  NU  \ 

2  3iV 


0  (4.25) 


Signal-to-Noise  plus  Interference  Ratio 

Far{^}  =  Var\r]]  +  Var\g^ 

=  yoZ:+yy5Zl 

4  trs  6iV 


(4.26) 


SNIR 


Far{^} 


M  K-\ 

■IZ^ 


SV/7: 


TW^  M  ^,  -1  p 

-'^0  I  y  y  Pk 
2Tf 


(4.27) 


The  footprint  of  the  geostationary  satellite  is  very  small  eompared  to  its  distanee 
to  the  earth.  Therefore,  all  the  reeeived  signal  powers,  and  7^  ean  be  eonsidered  equal. 
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(4.28) 


SNIR  = 


1 

N,  ^MjK-l) 
2TP^  3N 


P,T  =  E, 

£ 

— ^  =  SNR  {Signal  to  Noise  Ratio) 


SNIR  = 


1 

1  M{K-l) 
2 -SNR  3N 


(4.29) 


Equation  (3.13)  gives  the  bit  error  probability  of  a  BPSK  system.  We  need  to 
replaee  the  SNR  expression  in  that  formula  with  the  Signal-to-Noise  plus  Interferenee 
Ratio  found  in  Equation  (4.29).  The  bit  error  probability  for  BPSK  modulation  with 
spreading  is: 


r 

\ 

2-SNR  3N  j 

(4.30) 


where 

=M{K-\) 
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Bit  Error  Probability  with  Uncoded  BPSK  for  various  number  of  users 


0  2  4  6  8  10  12  14  16  18  20 

SNR  (dB) 

Figure  17.  Bit  Error  Probability  for  Uncoded  BPSK 


As  seen  in  Figure  17,  the  performance  is  not  very  good  without  coding,  and  it 
degrades  rapidly  as  the  total  number  of  users  increase. 


B,  PERFORMANCE  ANALYSIS  WITH  FORWARD  ERROR  CORRECTION 
CODING 

1.  Performance  Analysis  with  Convolutional  Encoding 

The  bit  error  probability  of  a  communication  system  with  convolutional  coding  is 
given  in  [16]  as: 

P,<t.  APM) 

Pj  :  The  code  distance  spectrum. 

Ej  (f^) :  The  probability  of  error  in  the  pair  wise  comparison  of  metrics. 
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P^^d)  is  defined  differently  for  Hard  and  Soft  Deeision  Deeoding. 


a.  Hard  Decision  Decoding 


P,(d): 


d  f 

z  p\^-pr 

H+i)  V  y 
2 


ydjl, 


d-i 


,d  odd 


,d  even 


•  «  1  V  w 

t=— +1  ^  ^ 

2 


( 

A 

(^snir-r)  =  q 

tl 

1  iVo  |M(^-1) 

1) 

1-Sm  3N  ) 

R :  Code  Rate 


b.  Soft  Decision  Decoding 

Pfd)  =  Q(^d-SNIR-R) 


R :  Code  Rate 


(4.31) 


(4.32) 


(4.33) 


(4.34) 
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Figure  18.  BER  for  50  Users  and  Hard  Decision  Convolutional  Decoding  with  Different 

Rates  and  Constraint  Lengths  {K^) 

Lower  bit  error  rates  can  be  achieved  by  higher  code  rates,  as  seen  in  Ligure  18. 
However,  higher  rates  decrease  the  information  bit  rate  dramatically.  Rate  1/2  performs 
almost  as  well  as  rate  1/120  and  it  does  not  affect  the  information  bit  rate  as  much. 
Therefore,  in  situations  where  extreme  reliability  is  required,  lower  rates  can  be  used,  but 
for  normal  usage,  rate  1/2  gives  sufficient  performance.  The  spectrum  distances 
algorithm  used  to  obtain  the  coefficients  gives  perfect  results  for  rates  1/4,  1//20, 
1/40,  1/60,  1/120  [22].  Because  of  these  perfect  spectral  distance  coefficients,  rate  1/20 
performs  better  than  rate  1/16.  Also,  higher  constraint  lengths,  K^,  improve  the 

performance,  but  this  makes  the  decoder  hardware  more  complex.  An  increase  in  the 

from  8  to  9  improves  the  performance  about  1  dB,  except  rate  1/120,  where  the 
improvement  is  less. 
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The  number  of  users  greatly  affects  the  performance  of  the  system.  As  seen  in 
Figure  19,  increasing  the  number  of  users  from  50  to  200  degrades  the  bit  error 
probability  substantially.  The  optimum  number  of  users  can  be  around  100  because  it 
gives  the  additional  channel  capacity  without  excessively  degrading  the  performance. 


Pjj  for  Rate  =  1/2  for  Hard  Decision  Decoding 


SNR 

Figure  19.  BER  for  Rate  1/2  Hard  Decision  Decoding  with  Different  Number  of  Users 

Figure  20  shows  the  difference  in  the  performances  of  hard  decision  and  soft 
decision  decoding  algorithms.  The  soft  decision  performs  better  than  the  hard  decision 
decoding;  however,  its  hardware  is  more  complex. 
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P.  for  50  Users,  K  =  8,  N  =  256 
b  ’  c 
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Figure  20.  Differences  Between  Hard  Decision  and  Soft  Decision  Decoding 
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From  Equation  (4.32) 


Pjj  for  Rate  =  47/63  for  Reed-Solomon  Decoding 


n 
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Figure  21 .  for  Reed-Solomon  Deeoding 

As  mentioned  in  Chapter  111,  Reed-Solomon  eodes  are  not  very  effieient  for 
normal  distribution  of  errors.  RS  eodes  are  more  effective  against  bursts  of  errors.  Figure 
21  shows  that  the  performance  of  RS  decoding  is  not  as  effective  as  convolutional  coding 
when  compared  to  Figure  19. 
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3. 


Performance  Analysis  with  Concatenated  Coding:  RS-RS  Codes 


n,  =  I"*'  - 1 ,  n^-k^=  2^j 

^2  =  2”"  - 1 ,  n2-k2=  2^2 


k 

Inner  Code  Rate,  r,  =  — 


Outer  Code  Rate, 


Coneatenated  Code  Rate  R  =  r,n  = 


k,k. 


V2 


Symbol  Error  Probability  at  output  of  inner  deeoder: 


p.  <  E 


j  +  h 


y=/,+i 


VJ) 


Pj{^-Ps) 


«l-7 


Ps=^-{^-pT 


p  =  Q 


R 


\  2-SNR 


,  ^M(A-l) 


3A 


Bit  error  probability  at  output  of  outer  deeoder  (i.e.  Coneatenated  bit  error 
probability) 


p  <«2+l  V  J  +  h 


2n 


Z 


2  j=t2+l 


J  ) 
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Figure  22.  Concatenated  RS-RS  Codes 

Concatenated  RS-RS  codes  do  not  perform  very  well  up  to  a  certain  SNR.  After 
that  point,  the  performance  increases  rapidly  with  slight  increases  in  SNR.  The  decrease 
in  the  performance  becomes  faster  as  the  number  of  users  increases. 


4.  Performance  Analysis  with  Concatenated  Coding:  Convolutional-RS 
Codes 

n,  =  2'”‘  - 1 ,  n^-k^=  2t^ 

^2  =  2”"  - 1 ,  ^2  “  ^2  ~ 
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Inner  Code  Rate, 


Outer  Code  Rate, 


n, 

K 


Coneatenated  Code  Rate  R  = 

The  bit  error  probability  at  the  output  of  the  inner  eonvolutional  deeoder; 


d^dfree 

The  performanee  of  only  soft  deeision  inner  deeoding  will  be  examined  in  this 

thesis. 


P^{d)  =  Q[dd-SNIR-R') 

The  symbol  error  probability  at  the  input  of  outer  RS  deeoder: 

/  \^2 

The  bit  error  probability  at  the  output  of  outer  deeoder  (Coneatenated  Bit  Error 
Probability): 


<zvtl  ^ 


^^2  j=t2+l  ^2 


J  J 


pA'-p..) 


n2-J 
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P|j  for  Concatenated  Convolutional  -  RS  codes 


Figure  23.  Concatenated  Convolutional  -  RS  codes 


Concatenated  systems  do  not  always  perform  better  than  the  single  coding 
structures.  When  Figure  23  is  compared  to  Figure  19,  we  can  clearly  see  that  rate  1/2 
convolutional  coding  performs  better  than  the  concatenated  system  which  has  an  inner 
code  rate  of  3/4  and  outer  code  rate  of  43/63.  The  overall  code  rate  is  0.5119,  which  is 
very  close  to  rate  1/2.  However,  the  concatenated  system  can  correct  both  normally 
distributed  errors  and  burst  errors.  The  concatenated  Convolutional  -  Reed-Solomon 
codes  are  affected  by  the  increase  in  the  number  of  users  worse  than  the  other  codes  as 
seen  in  Figure  24. 
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SNR 


Figure  24.  Pi^  for  Concatenated  Convolutional-RS  codes  with  different  number  of  users 

D.  CONCLUSION 

We  have  given  a  detailed  analysis  of  a  spread-spectrum  communication  system 
with  various  coding  schemes.  Each  analysis  has  different  parameters  that  can  be  changed 
in  order  to  find  the  optimum  satellite  communication  system.  From  the  results  of  this 
analysis,  the  suggested  VS  AT  system,  which  utilizes  CDMA,  will  use  concatenated 
convolutional-RS  coding,  support  100  users  and  still  provide  low  error  rates.  Actually,  a 
system  with  50  users  performs  better,  but  the  system  capacity  is  very  low  in  that  case. 
The  concatenated  system  does  not  perform  as  well  with  the  convolutional  coding  alone, 
but  concatenation  provides  extra  protection  for  burst  errors,  where  convolutional  coding 
is  effective  only  against  normally  distributed  errors.  Finally,  a  system  with  inner 
convolutional  code  rate  of  3/4  and  outer  RS  code  rate  of  57/63  gives  the  best  result  of 

tor  any  — ^  . 

N, 
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V.  CONCLUSIONS 


In  this  study,  a  satellite  communieation  system,  which  uses  VSATs  as  the  ground 
terminals,  was  examined.  The  system  utilizes  CDMA  with  Forward  Error  Correction.  The 
CDMA  system  uses  PN  sequences  and  Walsh  Codes  to  generate  a  spread  signal.  Each 
VSAT  uses  a  different  PN  sequence.  Each  user  in  a  VSAT  uses  different  Walsh  Codes. 
Since  Walsh  codes  are  orthogonal,  interference  among  the  users  is  eliminated.  PN 
sequences  are  not  completely  orthogonal;  therefore  interference  exists  due  to  other 
VSATs.  To  reduce  the  effects  of  interference  and  channel  errors,  the  spread  signal  is 
encoded  with  different  error  correction  mechanisms  and  its  performance  was  analyzed. 
The  results  for  bit  error  probabilities  are  upper  bounds,  not  the  exact  results.  Particularly, 
Convolutional  and  Reed-Solomon  codes  were  examined. 

Convolutional  codes  provide  low  bit  error  rates  and  they  are  flexible.  Good 
convolutional  codes  provide  the  maximum  coding  gain.  The  results  for  different  code 
rates  show  that  higher  code  rates  result  in  lower  bit  error  rates.  Rate  1/120  has  the  best 
performance  among  all,  but  this  rate  decreases  the  information  rate  dramatically.  On  the 
other  hand,  rate  1/2  performs  almost  as  good  as  rate  1/120  and  it  does  not  reduce  the 
information  rate  that  much.  In  cases  where  extremely  reliable  communications  are 
needed,  lowering  the  data  rate  and  increasing  the  code  rate  can  satisfy  the  users. 
However,  lowering  the  data  rate  makes  voice  communication  impossible  after  a  certain 
level.  In  these  cases,  only  data  communication  can  be  established. 

Reed-Solomon  (RS)  codes  are  used  in  many  modern  communication  systems. 
They  are  effective  for  bursts  of  errors,  and  they  are  not  very  good  in  correcting  the 
random  errors.  The  increase  in  the  error  rate  of  RS  codes  is  more  than  convolutional 
codes  for  the  same  amount  of  increase  in  the  number  of  users.  The  block  length  for  the 
RS  code  in  this  study  is  63  ( 2'’  - 1 ).  Eor  lower  bit  error  rates,  longer  block  lengths  can  be 
selected. 

The  optimum  results  are  obtained  by  using  concatenated  codes.  Concatenated 
codes  are  used  for  deep  space  missions  by  the  NASA.  Usually,  the  inner  codes  are 
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Viterbi  decoded  convolutional  codes  and  the  outer  codes  are  Reed-Solomon  block  codes. 
The  inner  code  rate  used  in  this  study  was  always  y4,  because  this  rate  has  a  good 
performance  and  it  does  not  decrease  the  information  rate  very  much.  The  outer  code  can 
be  changed  according  to  the  operational  requirements.  An  inner  code  rate  of  3/4  together 
with  an  outer  code  rate  code  rate  43/63  give  a  concatenated  code  rate  0.5119.  When  this 
result  is  compared  to  rate  1/2  convolutional  code,  the  bit  error  probability  of  the  former  is 
not  very  good.  For  higher  Signal-to-Noise  ratios,  concatenated  codes  have  a  greater 
advantage  over  convolutional  codes. 

The  VS  AT  system  that  can  be  proposed  based  on  this  study  is  a  DS-CDMA 
system  with  rate  3/4  Viterbi  decoded  inner  convolutional  code  and  rate  57/63  outer  RS 
code.  Using  an  interleaver  between  the  inner  and  outer  codes  can  further  improve  the 
burst  error  performance.  Also  using  different  code  rates  both  for  the  inner  and  the  outer 
codes  can  help  to  find  the  lowest  bit  error  rate  for  a  given  requirement. 


54 


LIST  OF  REFERENCES 


[1]  COL  Hook,  J.A.  Jr.,  Military  Dependence  on  Commercial  Satellite 
Communications  Systems  -  Strength  or  Vulnerability?  Air  War  College,  Air 
University,  Maxwell  Air  Foree  Base,  Alabama,  April  1999. 

[2]  Hadjitheodosiou,  M.H.,  F.P.  Coakley,  “Performanee  Optimization  of 
Multiple  Aecess  Protoeols  for  Multi-Serviee  VS  AT  Satellite  Networks,” 
IEEE  Telecommunications  Conference  Publication,  No.  404,  pp.  9-13, 
1995. 

[3]  Asgari  V.,  M.  Shiva,  “Channel  Coding  for  VSAT  Systems,”  IEEE 
Proceedings,  Vol.  142,  No.  5,  pp.  333-339,  Oetober  1995. 

[4]  Tzeng,  F.F.,  “Error  Protection  for  Satellite  Broadcast  of  HDTV  and 
Conventinal  TV,”  IEEE  Publications  0-7803-0917-0/93,  pp.  1617-1621, 
1993. 

[5]  Moheb,  H.,  C.  Robinson,  J.  Kijeski,  “Design  and  Development  of  Co- 
Polarized  Ku-band  Ground  Terminal  System  for  VSAT  Application,”  IEEE 
Publications  0-7803-5639-X/99,  pp.  2158-2161,  1999. 

[6]  Jonathon,  Y.,  C.  Cheah,  M.E.  Davis,  “Tone  Interference  in  VSAT 
Sptectrum,”  IEEE  Transactions  on  Communications,  Vol.  45,  No.  9,  pp. 
1035-1038,  September  1997. 

[7]  Rana,  Hamid  AA,  J.  McCoskey,  W.  Check,  “VSAT  Technology,  Trends, 
and  Applications,”  IEEE  Proceedings,  Vol.  78,  No.  3,  pp.  1087-1095,  July 
1990. 

[8]  Tamar  A.  Mehuron,  “Space  Almanac,”  Air  Force  Magazine,  Vol.  84,  No.  8, 
pp.  27-57,  August  2001. 

[9]  Evans,  B.G.,  Satellite  Communication  Systems,  3'^‘*  ed..  United  Kingdom. 
The  Institution  of  Electrical  Engineers,  1999. 


55 


[10]  Hughes,  C.D.,  C.  Soprano,  F.  Feliciani,  “Satellite  Systems  in  a  VSAT 
Environment,”  Electronics  &  Communication  Engineering  Journal,  Vol.  5, 
No.  5,  pp.  285-291,  October  1993. 

[11]  Maral.  G.,  M.  Basque!,  Satellite  Communications  Systems,  ed.,  Toulouse, 
France.  Wiley,  1998. 

[12]  Roddy,  D.,  Satellite  Communications,  3'^‘*  ed..  New  York.  McGraw  Hill, 

2001. 

[13]  Shannon,  C.E.,  “A  Mathematical  Theory  of  Communication,”  Bell  Systems 
Technology  Journal,  vol.  27,  pp.  623-656,  July  1948. 

[14]  Elbert,  B.R.,  The  Satellite  Communication  Ground  Segment  and  Earth 
Station  Handbook,  Artech  House,  2000. 

[15]  Rappaport,  T.S.,  Wireless  Communications,  2“‘*  ed.,  Prentice  Hall,  2002. 

[16]  Proakis,  J.G.,  Digital  Communications,  4**'  ed.,  McGraw  Hill,  2000. 

[17]  Ziemer,  R.E.,  R.E.  Peterson,  Introduction  to  Digital  Communications,  2"‘* 
ed.,  Prentice  Hall,  2000. 

[18]  Wicker,  S.B,  K.B.  Vijay,  Reed  Solomon  Codes  and  Their  Applications, 
IEEE  Press,  1994. 

[19]  Tighe,  J.E.,  Modeling  And  Analysis  of  Cellular  CDMA  Forward  Channel, 
Ph.D.  Dissertation,  Naval  Post  Graduate  School,  Monterey,  California, 
March  2001. 

[20]  Eee,  J.S.,  E.E.  Miller,  CDMA  Systems  Engineering  Handbook,  Artech 
House,  1998. 

[21]  ETJG.  Aras,  E.,  Title  of  Aras’  Thesis,  Naval  Post  Graduate  School, 
Monterey,  California,  September  2002 

[22]  Erenger,  P.,  P.  Orten,  T.  Ottoson,  “Code-Spread  CDMA  Using  Maximum 
Eree  Distance  Eow-Rate  Convolutional  Codes”,  IEEE  Transactions  on 
Communications,  Vol.  48,  No.  1,  pp.  135-144,  January  2000. 


56 


[23]  Clark,  G.C.Jr.,  J.B.  Cain,  Error  Correction  Coding  for  Digital 
Communications,  Plenum  Press,  New  York,  1981. 


57 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


58 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Teehnieal  Information  Center 

Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 

Naval  Postgraduate  Sehool 
Monterey,  California 

3.  Professor  Tri  T.  HA 

Naval  Postgraduate  School 
Monterey,  California 

4.  Professor  Herschel  H.  Loomis,  Jr 

Naval  Postgraduate  School 
Monterey,  California 

5.  Professor  Dan  C.  Boger 

Naval  Postgraduate  School 
Monterey,  California 

6.  Kara  Kuvvetleri  Komutanligi 

Kutuphane 

Bakanliklar,  Ankara,  TURKEY 

7.  Kara  Harp  Okulu  Komutanligi 

Kutuphane 

Bakanliklar,  Ankara,  TURKEY 

8.  Deniz  Harp  Okulu  Komutanligi 

Kutuphane 

Tuzla,  Istanbul,  TURKEY 

9.  Hava  Harp  Okulu  Komutanligi 

Kutuphane 

Tuzla,  Istanbul,  TURKEY 


59 


